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Physics-Based  Modeling  for  WUI  Fire  Spread  - Simplified 
Model  Algorithm  for  Ignition  of  Structures  by  Burning 
Vegetation 


David  D.  Evans,  Ronald  G.  Rehm,  and  Elisa  S.  Baker 
Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 


1.0  ABSTRACT 

A simple  physics-based  mathematical  algorithm  is  described  for  estimating  the  separation 
distances  necessary  to  prevent  radiative  ignition  of  structures  by  burning  vegetation.  The 
algorithm  is  directed  at  use  in  applications  such  as  the  USDA  Forest  Service  EcoSmart 
software  suite  of  tools.  The  algorithm  represents  a first  step  in  quantifying  the  fire  spread 
phenomena  involving  both  structures  and  vegetation  typical  of  wildland-urban  interface 
(WUI)  areas  and  residential  housing  communities.  Results  of  initial  experiments  to 
characterize  the  heat  release  rate  (HRR),  flame  height,  and  duration  of  burning  for 
individual  Douglas-fir  trees  are  provided.  In  addition,  research  required  to  understand 
better  and  quantify  wildland-urban  and  residential  community  fires  are  discussed. 


2.0  INTRODUCTION 

The  Center  for  Urban  Forest  Research,  the  Pacific  Southwest  Research  Center  of  the  U.S. 
Department  of  Agriculture  (USDA)  Forest  Service,  together  with  other  partners,  is 
developing  a web-based  suite  of  software  tools,  called  EcoSmart.  EcoSmart  is  being 
developed  for  use  by  individual  homeowners  to  provide  guidance  in  assessing  the 
contributions  of  landscape  vegetation  to  home  energy  conservation,  land  hydrology,  and 
fire  spread.  EcoSmart  provides  means  to  quantify  the  advantages  and  disadvantages  of 
landscape  vegetation  relative  to  the  costs. 

As  a partner  in  this  effort,  with  partial  support  from  the  Pacific  Southwest  Research 
Station,  the  Fire  Research  Division  of  the  Building  and  Fire  Research  Laboratory  of  the 
National  Institute  of  Standards  and  Technology  (NIST)  developed  a simple  physics-based 
mathematical  model  for  estimating  whether  a structure  will  be  ignited  by  nearby  tree 
crown  fires.  This  model  quantifies  the  burning  of  vegetation  in  a manner  that  has  not 
been  done  before.  The  model  developed  treats  each  tree,  shrub,  and  house  as  an 
individual  fuel  element.  In  order  to  predict  fire  spread  from  burning  fuel  elements  by 
ignition  of  other  fuel  elements,  the  ignition  and  burning  characteristics  of  each  fuel 
element  and  the  interactions  between  elements  has  to  be  defined  or  determined.  This 
creates  a need  for  ignition  and  burning  data  that  has  not  been  generally  measured  before 
for  vegetation  and  structures.  With  this  data  and  appropriate  models,  the  importance  of 
individual  fuel  elements  in  fire  propagation  can  be  assessed.  Resolving  individual  fuel 
elements  in  a model  is  necessary  in  order  to  properly  evaluate  the  WUI  and  residential 


housing  community  fires.  In  these  areas,  structures  are  an  important  asset  to  be  protected 
from  fire,  but  can  be  an  important  source  of  fire  spread  if  ignited.  Fire  spread  in  the  WUI 
or  residential  housing  communities  cannot  be  characterized  using  the  methods  based  on 
fuel  bed  classifications  that  are  used  widely  to  analyze  the  movement  of  fire  fronts  in 
forest  fires.  Regarding  each  structure  as  a fuel  element,  its  likely  fate  in  a fire  can  only  be 
determined  by  examining  its  possible  interactions  with  surrounding  fuel  elements  that 
can,  if  burning,  threaten  the  structure.  The  algorithm  developed  in  this  study  is  the  first 
step  in  creating  the  foundation  for  a new  method  to  analyze  fire  threat  in  WUI  areas  and 
residential  communities. 

Version  1.0  of  the  fire  algorithm  developed  by  NIST  is  available  for  and  compatible  with 
the  level  of  detail  in  other  portions  of  the  EcoSmart  software  under  development.  The 
algorithm  documented  in  this  report  represents  a first  step  in  quantifying  fire  spread 
phenomena  involving  both  structures  and  vegetation  typical  of  wildland-urban  areas  and 
residential  community  properties.  The  report  provides  new  data  on  HRR,  flame  height, 
and  duration  of  burning  for  individual  Douglas-fir  trees.  In  addition,  it  also  contains  a 
discussion  of  research  necessary  to  better  understand  and  quantify  wildland-urban  and 
residential  community  fire  spread. 


3.0  BACKGROUND 

The  development  of  the  EcoSmart  suite  of  software  tools  began  at  the  U.S.  Department  of 
Agriculture  Forest  Service,  Pacific  Southwest  Research  Station  in  2002.  There  was  a 
need  to  develop  an  algorithm  that  accounted  for  the  increased  risk  for  home  destruction 
by  fire  due  to  the  placement  of  trees  and  other  ignitable  vegetation  close  to  the  structure. 
This  factor  was  in  opposition  to  the  benefits  these  plantings  may  have  for  energy 
conservation  and  soil  hydrology.  Development  of  the  EcoSmart  software  suite  of  tools 
was  undertaken  to  allow  property  owners  to  evaluate  advantages  and  disadvantages, 
including  fire  safety  and  economics,  of  planting  landscaping  vegetation  around  homes. 

NIST  was  asked  by  the  Pacific  Southwest  Research  Station  to  develop  a physics-based 
algorithm  for  ignition  of  structures  by  landscape  vegetation  that  could  be  programmed 
into  the  EcoSmart  suite  of  software  tools  and  would  be  consistent  with  level  of 
information  and  desired  speed  of  operation  of  all  the  other  parts  of  EcoSmart.  Very  little 
information  was  found  in  the  literature  that  could  be  used  to  construct  a physics-based 
model  for  ignition  of  structures  by  burning  vegetation  of  the  type  and  arrangement  that 
would  be  present  on  landscaped  properties.  Principally,  there  was  a lack  of  data  on 
burning  characteristics  for  individual  trees  and  shrubs. 

Fire  is  an  extremely  complex  natural  phenomenon  [1-10].  Decades  of  research  to 
understand  and  quantify  fire,  together  with  revolutionary  advances  in  computer 
technology  and  computational  algorithms,  have  lead  to  the  development  of  sophisticated 
mathematical  models  that  can  predict  and  visualize  growth  and  spread  of  fire  under  a 
variety  of  conditions  [11-16].  However,  these  models  generally  require  substantial 
computational  and  data  resources.  Furthermore,  since  they  have  not  been  oriented  toward 
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the  prediction  of  WUI  or  residential  community  fires,  they  will  not  soon  be  available  for 
application  by  end  users  without  a substantial  research  investment  over  the  next  few 
years. 


3.1  PREVIOUS  MODELING  OF  WILDLAND  FUELS  AND  FIRES 

Fires  in  the  western  U.S.  have  been  headlines  in  the  news  for  several  recent  years. 
Changes  in  the  management  of  forested  lands,  the  increasing  intrusion  of  man  into  more 
remote  areas  and  the  cyclic  dry  periods  produced  by  El  Nino  - La  Nina  have  all  increased 
the  loss  of  man-made  structures  to  fires.  Wildland  fires  now  often  spread  into  the  built 
environment  causing  injury,  death  and  property  damage.  Concern  about  wildland  fires 
and  their  effects  on  man  and  the  built  environment  has  generated  increased  interest 
recently,  e.g.,  Livingston  [17],  Keller  [18],  Platt  [19].  While  the  need  to  address  fire 
spread  in  a mixed  environment,  containing  both  structures  and  wildland  fuels,  has  been 
acknowledged,  most  research  has  continued  to  focus  on  wildland  fires.  Operational 
models  that  predict  fire  spread  in  wildland  fuels  are  well  developed  and  widely  used,  but 
do  not  yet  address  WUI  fires. 

For  over  thirty  years,  modeling  of  fire  dynamics  in  wildland  fuels  has  been  an  ongoing 
research  and  development  activity  [20-26].  It  has  lead  to  practical  or  “operational” 
models  that  are  regularly  used  by  land  managers  to  predict  fire  growth  and  spread  for 
both  planned  or  “controlled”  bums  and  unexpected  wildland  fires.  These  operational 
models  have  served  land  managers  and  firefighters  well  both  for  planning  controlled 
bums  and  suppression  of  wildfires.  In  the  United  States,  the  mathematical  models 
BEFIAVE  [11]  and  FARSITE  [12]  are  two  of  the  most  widely  used  operational  models. 

Both  are  based  upon  the  pioneering  research  of  Rothermal  [20,  21],  who  developed  the 
framework  for  these  models  and  determined  many  of  the  empirical  relations  used  by  it.  In 
these  models,  wildland  fuels  are  characterized  by  13  standard  fire  behavior  fuel  types. 
These  fuel  characterizations  and  their  behavior  in  w ildland  fires  is  the  subject  of  several 
of  the  books  cited  above:  Chandler  et  al  [1],  Chandler  et  al  [2],  Luke  and  McArthur  [3], 
Pyne  et  al  [4],  and  Brown  and  Davis  [5]  lists  the  13  types  and  their  characteristics. 

The  basic  model  regards  all  of  the  fuel  to  cover  the  land  surface  as  though  it  were 
"painted”  on;  i.e.,  all  of  the  combustible  material  resides  in  fuel  beds  that  have  no  vertical 
extent,  but  are  strictly  two  dimensional.  Fire  then  spreads  over  these  fuel  beds  as  a 
contour  which  divides  unbumed  fuel  from  burned  fuel  or  combustion  products.  Fuel  beds 
can  be  homogeneous  or  inhomogeneous,  isotropic  or  anisotropic,  and  they  can  be 
composed  of  a single  fuel  type  or  a combination  of  fuel  types.  Furthermore,  the  land 
surface  may  be  flat  or  may  have  significant  topographical  features,  and  all  of  these 
variations  make  for  differences  in  how  the  fire  line  progresses. 

Inclusion  of  improved  descriptions  of  the  physics  in  these  models  has  progressed  slowly 
in  recent  years  and  has  not  taken  advantage  of  the  revolution  in  computer  hardware,  and 
computational  algorithms  and  software  that  has  occurred  over  this  same  time  frame. 


By  comparison,  over  the  past  few  decades,  many  of  the  physical  and  chemical  processes 
occurring  in  fires  have  become  much  better  understood,  and  physics-based  mathematical 
and  computational  models  have  been  developed  to  take  advantage  of  this  improved 
understanding.  One  such  example  is  the  work  of  Clark  [27]  and  colleagues  in  which  the 
coupling  between  large-scale  wildland  fires  and  the  atmosphere  has  been  studied. 

Another  is  modeling  of  the  growth  and  spread  of  fires  in  structures,  which  has  advanced 
dramatically  using  the  methods  of  computational  fluid  dynamics  (CFD)  and,  more 
generally,  mathematical  and  computational  methods  that  may  be  described  as 
computational  mechanics  and  computational  combustion  [28-31]. 

For  models  of  combustion,  the  usual  length  scales  of  interest  are  measured  in  very  small 
fractions  of  a meter  (millimeters  or  less)  to  meters.  Alternately,  for  problems  related  to 
computational  mechanics  or  CFD,  the  length  scales  of  interest  can  vary  between  meters 
and  tens  or  even  hundreds  of  meters.  These  scales  should  be  contrasted  with  models  for 
wildland  fires,  where  the  lengths  are  measured  in  kilometers  or  even  hundreds  of 
kilometers. 

3.2  PREVIOUS  MODELING  OF  WUI  FIRES 

Very  little  progress  has  been  made  on  understanding  and  predicting  fire  spread  at  the 
Wildland-Urban  Interface  (WUI),  where  both  wildland  fuels  and  man-made  structures 
coexist.  At  the  WUI,  the  length  scales  are  intermediate  to  the  wildland  scales  on  the  one 
hand  and  to  those  associated  with  individual  structures  on  the  other.  As  a result,  a model 
for  a WUI  fire  must  describe  individual  structures  and  trees  in  some  detail,  and  must  be 
able  to  determine  fire  spread  between  these  individual  fuel  sources,  a much  more 
challenging  task!  One  of  the  few  studies  to  attack  this  very  difficult  problem  is  the  work 
of  Maranghides  [32],  who  concluded  that  there  was  an  extreme  lack  of  relevant  data. 

J.  D.  Cohen  [33,  34]  advocates  the  concept  of  a Home  Ignition  Zone,  defined  as  the  home 
and  its  surrounding  area  where  heat  transfer  by  radiation  and/or  convection  from  burning 
vegetation  can  ignite  it.  He  is  developing  a model  called  the  Structure  Ignition 
Assessment  Model  (SIAM)  to  quantify  this  concept,  which  is  being  recognized  as 
relevant  to  homeowners  and  firefighters.  The  model  estimates  conditions  for  home 
ignition.  Once  a home  ignites  in  the  calculation,  it  is  considered  lost.  SIAM  does  not 
address  the  contribution  of  the  burning  structure  to  subsequent  fire  spread  in  the  area.  It 
therefore  has  limited  application  in  the  WUI,  especially  in  cases  of  housing 
developments,  where  homes  may  be  located  within  the  Home  Ignition  Zone  of  its 
neighbors. 


3.3  NIST  RESEARCH  CONTRIBUTION 

In  order  to  understand  better  the  burning  of  landscape  trees,  a series  of  burns  was 
performed  in  the  NIST  large  fire  facility  to  measure  the  burning  characteristics  of 
individual  Douglas-fir  trees.  The  measured  burning  characteristics  of  the  trees  were  used 
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as  a basis  for  the  construction  of  the  algorithm  for  ignition  of  structures  by  burning 
vegetation. 

During  this  research,  three  papers  were  prepared  that  provide  background  and  discussion 
of  fire  spread  in  the  wildland-urban  interface  and  in  residential  housing  communities  [35- 
37].  In  recent  WUI  fire  disasters,  visual  evidence  suggests  that  fully  involved  burning 
houses  contribute  to  fire  spread  and  destruction.  Examples  are  given  to  show  that  the 
long-duration  and  high-intensity  burning  of  a structure  can  contribute  to  altering  the  fire- 
generated winds  [36].  In  addition,  early  implementations  of  the  structure  ignition 
algorithm  in  EcoSmart  are  illustrated  [35]. 

4.0  BURNING  CHARACTERISTICS  OF  DOUGLAS-FIR  TREES 

Although  landscaping  trees  and  shrubs,  if  ignited,  can  endanger  near-by  structures,  very 
little  quantitative  information  exists  about  ignition  and  burning  characteristics  of 
landscape  vegetation.  Generally,  experiments  carried  out  to  measure  the  HRR  of  trees 
have  been  conducted  to  assess  the  hazards  of  Christmas  tree  fires  [38-40],  No 
information  has  been  found  to  quantify  how  the  HRR,  bum  duration,  and  flame  height 
vary  with  the  size  of  a tree  or  tree  species.  Two  studies  were  found  in  which  the  burning 
characteristics  of  shrubs  were  measured  [41,  42], 

Previously,  the  HRR  of  dry  Scotch  Pine  trees  was  measured  using  the  instrumented  hood 
in  the  Large  Fire  Laboratory  at  NIST  [38],  However,  the  study  included  only  one  size 
tree.  The  same  is  true  for  the  study  of  Douglas-fir  tree  burning  carried  out  by  Babrauskas 
[39].  To  take  advantage  of  existing  data,  in  this  study  we  continued  to  bum  cut  Douglas- 
fir  trees.  Three  different  commercial  sizes  were  selected  for  the  experiments.  The  height 
of  the  trees  was  limited  by  the  HRR  capacity  of  the  Large  Fire  Laboratory  exhaust  hood. 
The  tree  bums  provided  important  data  for  the  algorithm  development.  Further 
experiments  are  needed  to  fully  understand  the  burning  characteristics.  Data  obtained 
from  the  bums  provided  information  for  the  development  of  the  algorithm  that  was  not 
available  from  any  other  source. 

4.1  EXPERIMENTS 

Prior  to  burning,  the  cut  trees  were  allowed  to  dry  for  17  days  as  discussed  below.  All 
tree  bums  were  conducted  during  a single  week  in  the  Large  Fire  Laboratory.  This 
included  time  needed  for  instrument  set-up  and  mounting  each  tree  for  burning.  During 
each  bum,  the  HRR,  mass  loss,  flame  height,  vertical  temperature  profiles  and  radiant 
heat  flux  emitted  from  the  fire  were  measured. 

4.1.1  Tree  Preparation 

The  Douglas-fir  trees  were  purchased  from  a local  nursery  that  grows  Christmas  trees. 
Three  trees  each  with  commercial  sale  heights  of  1.2  m,  2.4  m,  and  3.7  m (4  ft,  8 ft,  and 
12  ft)  were  cut  on  July  6,  2003.  These  trees  were  delivered  to  NIST  the  following  day 
and  weighed.  As  the  trees  are  a large  fuel  load,  for  fire  safety  considerations,  the  trees 
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were  dried  on  racks  outside  the  Large  Fire  Laboratory.  The  racks  were  covered  by  tarps 
to  shield  the  trees  from  rain.  After  10  days,  the  trees  were  brought  inside  and  weighed 
again  to  monitor  the  moisture  lost  to  evaporation.  Some  of  the  trees  had  the  base  of  the 
trunks  trimmed  to  create  a flat  base  for  mounting  on  a stand  that  supported  the  tree 
upright  for  burning. 

The  trees  were  photographed  and  heights  were  measured  (see  Table  1 ).  The  crown 
height  is  the  height  from  the  bole  (height  of  the  first  live  branches)  to  the  top  branch  as 
shown  in  Figure  1 . This  differs  from  the  commercial  height  of  the  tree  which  is 
measured  from  the  base  to  the  tip  of  the  stem.  The  crown  height  best  represents  the 
burnable  fuel  height  by  eliminating  the  top  of  the  stem  and  base  of  the  tree  up  to  the  bole 
height.  The  stem  needle  growth  on  the  test  Douglas-fir  trees  that  was  eliminated  from  the 
burnable  mass  by  this  definition  extend  0.3  m above  the  highest  branch  in  some  cases, 
but  only  added  a small  amount  to  the  burnable  fuel  mass. 
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Tree 

Test  No. 

Crown 

Height 

(m) 

Bole 

Height 

(m) 

Total 

Height 

(m) 

Weight 

7/7/03 

(kg) 

Weight 

7/17/03 

(kg) 

Weight 

Pre-burn 

(kg) 

Weight 

Pre-burn 

Initial 

(7/7/03) 

4-1 

Test  1 

1.37 

0.08 

1.45 

5.16 

3.014 

2.85 

0.55 

4-2 

Test  4 

1.30 

0.19 

1.49 

3.54 

1.607 

1.549 

0.44 

4-3 

Test  7 

1.42 

0.22 

1.64 

8.4 

4.515 

4.372 

0.52 

8-1 

Test  2 

2.31 

0.27 

2.58 

23.11 

17.597 

17.597 

0.76 

8-2 

Test  5 

2.62 

0.43 

3.05 

24.4 

18.103 

17.637 

0.72 

8-3 

Test  8 

2.44 

0.48 

2.92 

40.83 

29.61 

28.751 

0.70 

12-1 

Test  3 

3.10 

0.69 

3.78 

38.17 

30.116 

29.92 

0.78 

12-2 

Test  9 

3.20 

0.71 

3.91 

41.26 

33.08 

31.92 

0.77 

12-3 

Test  6 

3.33 

0.50 

3.82 

38.34 

29.345 

29.14 

0.76 

Table  1 Douglas-fir  tree  identifications  and  physical  measurements 


4.1.2  Large  Fire  Laboratory 

The  trees  were  burned  in  the  NIST  Large  Fire  Laboratory  located  on  the  NIST  site  in 
Gaithersburg,  Maryland.  The  facility  has  several  instrumented  hoods  that  employ  oxygen 
consumption  calorimetry  to  measure  the  HRR.  The  largest  calorimeter  hood,  which  was 
used  for  the  tree  bums,  has  a design  capacity  for  HRRs  of  10000  kW.  The  hood 
measures  9 m x 12  m and  the  lower  edge  of  the  hood  is  4.5  m above  the  floor  level.  At 
its  highest,  in  the  center  of  the  hood,  the  clear  space  above  the  floor  is  7.9  m.  The 
products  of  combustion  are  collected  by  the  hood  and  flow  though  ducts  for  measurement 
and  then  are  exhausted  from  the  test  building.  After  leaving  the  building  the  fire  products 
are  treated  by  a pollution  abatement  system  before  being  discharged  into  the  atmosphere. 

4.1.3  Instrumentation 

The  experiments  were  video  taped  from  three  different  perspectives,  and  still  digital 
photographs  were  also  taken.  The  tree  and  stand  were  on  a platform  supported  by  a load 
cell  that  measured  the  mass  loss  as  a function  of  time.  Additionally,  a much  larger 
platform  underneath  the  first  was  supported  by  load  cells  at  each  of  the  four  comers,  to 
verify  the  accuracy  of  the  first  load  cell,  and  serve  as  a backup. 

Additional  instrumentation  included  the  following: 

Seven  thermocouples  were  mounted  on  a 4.3  m (14  ft)  vertical  thermocouple  support 
stand  placed  approximately  a meter  from  each  test  tree.  These  thermocouples  were 
mounted  on  horizontal  arms  at  0.61  m (2  ft)  intervals,  so  that  the  sensing  tips  would  be 
near  the  centerline  of  the  tree  crown  in  the  buoyant  gas  nearest  to  the  stem.  The 
thermocouples  were  3 mm  diameter  stainless  steel  sheathed,  Type  K Chomel-Alumel, 

0.91  m (3  ft)  long.  Voltage  output  from  the  thermocouple  array  was  recorded  on  the 
Large  Fire  Laboratory  data  acquisition  system. 

The  first  three  tests  had  4 water-cooled  total-heat-flux  gauges,  which  were  mounted  on  a 
second  4.3  m (14  ft)  vertical  support  located  3 m (9.8  ft)  from  the  stem  of  the  tree.  The 
first  gauge  was  located  0.61  m (2  ft)  above  the  ground  level,  and  the  remaining  three 
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gauges  were  spaced  vertically  every  1 .22  m (4  ft).  Additionally,  two  heat  flux  gauges 
were  placed  aiming  upward  at  ground  level  1 .5  m (4.9  ft)  and  3 m (9.8  ft)  from  the  tree 
stem.  For  the  fourth  through  ninth  tests,  a seventh  heat  flux  gauge  was  added,  10  m (32.8 
ft)  away  from  tree  stem.  This  gauge  was  angled  at  20°  to  view  the  entire  flame. 

4.2  MEASUREMENTS 

Figure  2 presents  a timeline  of  images  from  the  video  record  of  all  Douglas-fir  tree  burns 
except  Test  9 (tree  12-2)  that  was  lost.  Each  tree  was  ignited  at  the  base  of  the  crown 
using  a propane  torch  with  two  flames  for  simultaneous  ignition  of  the  crown  on  opposite 
sides  of  the  stem.  Depending  on  the  moisture  content,  some  trees  required  5 s of 
torch  application  to  ignite,  others  required  15  s.  Once  ignited,  each  tree  was  consumed 
by  rapid  burning  in  approximately  60  s. 

All  trees  except  Tree  12-1  behaved  similarly  in  the  burn  tests.  Tree  12-1  (Test  3)  did  not 
get  fully  involved,  and  required  a re-ignition  to  sustain  burning.  Because  of  this  unusual 
burning  behavior  the  results  for  Test  3 are  not  included  in  the  analysis.  As  there  were 
two  other  trees  of  the  same  nominal  size  remaining  in  the  data  set,  eliminating  Test  3 
from  the  analysis  was  regarded  as  not  significant. 
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Tree  4-1 


Figure  2a.  Pictures  every  5 s of  the  Douglas-Fir-Tree  Burns  starting  from  the  application 
of  the  igniter  until  55  s after  that. 
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Figure  2b.  Pictures  of  the  Douglas-Fir-Tree  Burns  for  times  between  60  s and  110  s after 
application  of  the  igniter. 
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4.2.1  Flame  Heights 


The  mean  flame  height,  Htiame  was  determined  for  each  tree  burn  by  analyzing  video  and 
still  photos  showing  both  the  burning  tree  and  the  installed  thermocouple  supports.  The 
flame  height  (see  figure  3)  was  chosen  to  be  the  distance  from  the  base  of  the  flame  near 
the  tree  bole  (the  lowest  branch  on  the  tree)  to  the  top  of  the  flame  w here  the 
intermittency  is  0.5.  A flame  intermittency  of  0.5  means  the  visible  flame  is  seen  half  of 
the  time.  For  each  bum,  still  photographs  were  used  to  determine  the  flame  height, 
guided  by  the  video  record  showing  flame  intermittency.  The  crown  height,  flame  height 
and  their  ratio  for  each  burn  are  shown  in  Table  2,  and  their  use  in  the  tree-bum  model  is 
presented  in  Sections  4.3.1  and  5.1.2. 


Tree 

Hcrown  (m) 

H flame  (m) 

Hflame  /Hcrown 

4-1 

1.37 

4.2 

3.06 

4-2 

1.30 

3.47 

2.68 

4-3 

1.42 

4.36 

3.07 

8-1 

2.31 

4.61 

1.99 

8-2 

2.62 

5.05 

1.93 

8-3 

2.44 

5.61 

2.30 

12-1 

3.10 

N/A 

N/A 

12-2 

3.20 

5.08 

1.59 

12-3 

3.33 

4.69 

1.41 

Table  2 Flame  heights 


4.2.2  Heat  Release  Rate  and  Mass  Loss 


Figure  4 shows  a plot  of  the  measured  HRR  for  each  of  the  tree-burns.  The  peak  HRR 
occurs  at  different  times  in  part  because  of  variation  in  ignition  times.  Peak  HRR  and 
weight  loss  measured  during  the  test  are  recorded  in  Table  3.  The  weight  loss  data  in 
Table  3 are  from  the  measurements  made  by  the  top  load  cell  located  directly  under  the 
tree. 
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Figure  4 HRR  as  a function  of  time 


Tree 

HRRpeak 

(kW) 

Tree  Weight 
Pre-burn 

(kg) 

Weight  Loss 
(kg) 

Weight  loss  / 
pre-bum  weight 

4-1 

780 

2.85 

1.44 

0.51 

4-2 

720 

1.55 

1.04 

0.67 

4-3 

1660 

4.37 

2.67 

0.61 

8-1 

3350 

17.60 

9.13 

0.51 

8-2 

3400 

17.64 

8.6 

0.49 

8-3 

5040 

28.75 

12.53 

0.44 

12-2 

5000 

31.92 

13.49 

0.42 

12-3 

4000 

29.14 

7.91 

0.27 

Table  3 Measured  peak  HRR  and  burned  weight  loss  for  tree 
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4.2.2  Other  Measurements 


As  part  of  the  large  scale  experiments,  many  measurements,  not  needed  immediately  for 
the  formulation  of  the  algorithm,  were  recorded.  These  included  temperature 
measurements  from  an  array  of  thermocouples  within  and  above  the  crown  of  the  tree, 
and  radiant  and  total  heat  flux  from  the  fire  at  various  distances  and  heights.  These  data 
will  be  reduced  and  analyzed  when  additional  resources  become  available. 

4.3  ANALYSIS 

The  series  of  tree-burn  experiments  provided  initial  data  to  start  the  formulation  of  the 
fire  source  portion  of  the  algorithm  to  predict  structural  ignition  from  burning  vegetation. 
Other  observations  from  the  data,  not  used  directly  in  the  development  of  the  algorithm, 
are  recorded  in  this  section  for  possible  use  in  other  models  for  wildland-urban  interface 
fires. 

4.3.1  Flame  Height 

The  mean  flame  heights  (Htiame)  were  compared  to  the  crown  height,  Hcr0Wn,  to  determine 
if  there  was  a correlation  between  a tree’s  height,  and  the  resulting  flame  height.  The 
results  are  shown  in  figure  5.  Generally  the  larger  trees  produced  a flame  height  tw  ice 
the  crown  height.  The  smallest  trees  produced  a flame  three  times  the  crown  height.  This 
data  is  also  influenced  by  differences  in  initial  moisture  content.  From  the  data  in  Table 
1,  the  smallest  trees  lost  about  one-half  of  the  initial  w eight  before  burning.  The  larger 
trees  lost  only  one-fourth  of  the  initial  weight  before  burning.  More  experiments  are 
needed  to  understand  the  effect  of  moisture  and  crow  n height  on  flame  height.  At  this 
time,  the  best  estimate  of  flame  height  for  application  to  landscape  trees  is  the  data  for  the 
larger  trees  with  greater  moisture.  That  is  Htiame  / HcroWn  = 2.0 


Figure  5 Flame  height  data 
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Tree 

H crown  (m) 

Hfiame  (m) 

H flamed  Hcrown 

4-1 

1.37 

4.2 

3.06 

4-2 

1.30 

3.47 

2.68 

4-3 

1.42 

4.36 

3.07 

8-1 

2.31 

4.61 

1.99 

8-2 

2.62 

5.05 

1.93 

8-3 

2.44 

5.61 

2.30 

12-2 

3.20 

5.08 

1.59 

12-3 

3.33 

4.69 

1.41 

Table  4 Tree  height  and  flame  height  data 


4.3.2  Heat  Release  Rate 

NIST  tree  burn  experiments  demonstrated  that  there  is  a similarity  in  the  growth  curves  to 
peak  HRR  amongst  all  the  trees,  regardless  of  size.  The  heat  release  data  from  the 
experiments  shows  wide  variation  in  the  time  to  peak  HRR  (Figure  3).  The  larger  trees 
took  longer  to  ignite,  and  have  a steeper  rate  of  increase,  but  once  the  fire  was  steadily 
increasing,  peak  HRR  was  obtained  in  approximately  30  s for  all  trees.  Figure  6 shows 
the  HRR  curves  aligned  such  that  the  peak  values  occur  at  30  s.  From  this  figure,  it  is 
apparent  that  the  size  of  the  tree  plays  a large  role  in  determining  the  peak  HRR,  but  has  a 
relatively  minor  affect  on  the  length  of  the  bum. 

Using  60  s as  a common  time  interval  for  burning,  the  HRR  of  the  burning  Douglas-fir 
trees  can  be  represented  well  by  a simple  isosceles  triangle  shape.  In  this  case  the  base  of 
the  triangle  is  60  s and  the  height  is  the  peak  HRR.  Figure  7 


Figure  6 Heat  release  rate  curves  aligned  at  their  peak  values. 


shows  this  approximation  as  it  applies  to  two  of  the  measured  HRR  curves.  This  result  is 
not  used  in  the  algorithm  developed  in  this  study,  but  may  be  an  important  and  useful 
approximation  in  future  modeling  of  crown  fires. 
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Figure  7 Triangle  approximation  for  the  heat  release  rate  of  Douglas-fir  crown  fires. 


The  burning  of  Douglas-fir  trees  has  also  been  studied  to  evaluate  the  hazards  of 
Christmas  trees  in  homes.  In  an  extensive  study  of  Christmas-tree  bums  by 
Babrauskas  [39],  peak  FIRR  were  correlated  with  total  mass  of  the  tree  and  the  moisture 
content  of  the  needles.  A plot  of  the  correlation  and  the  curve-fit  from  the  published 
results  of  that  study  are  reproduced  here  as  Figure  8 and  Equation  1 . 


Q peak  /maSS  = ^ 


5.84-0.01 7M 


(1) 


Where, 

q k is  the  peak  heat  release  rate  [kW], 

mass  is  the  total  pre-bum  mass  of  the  tree  [kg], 

M is  the  moisture  of  the  needles  [mass  fraction  or  percent  by  weight]. 

Because  these  results  contain  a quantification  of  the  effect  of  moisture  on  the  peak  HRR 
and  the  trees  burned  were  similar  to  those  burned  in  this  study,  the  results  published  by 
Babrauskas  will  be  used  in  the  formulation  of  the  algorithm.  As  Babrauakas  normalized 
his  peak  HRR  using  the  entire  mass  of  the  cut  Douglas-fir  tree,  his  formulation  will  need 
to  be  modified  to  w ork  with  the  burnable  mass  of  the  tree  estimated  from  experiments 
conducted  in  this  study. 
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Figure  8 - Peak  heat  release  rate  from  Douglas-fir  trees  as  a function  of  moisture 
and  mass  [39]  (reproduced  with  permission  from  the  Society  of  Fire 
Protection  Engineers) 

5.0  IGNITION  OF  STRUCTURES  ALGORITHM  FOR  ECOSMART 

In  wildland-urban  and  residential  community  areas,  fire  can  spread  rapidly,  often  driven 
by  wind,  to  threaten  homes  and  other  structures.  There  are  several  modes  by  which 
burning  vegetation  can  threaten  a home  [1-8].  For  example,  flames  from  burning  trees 
and  shrubs  close  to  a home  can  ignite  combustible  exteriors  by  direct  contact.  In  addition, 
small  fires  in  collected  debris,  such  as  dry  leaves  and  needles  dropped  on  and  around  the 
home,  can  start  fires  at  corners  and  in  small  crevices  of  the  structure.  Furthermore, 
ground  fires,  which  propagate  in  dead  debris  such  as  leaves  or  pine  needles  or  in  dry  live 
fuels  such  as  grasses,  can  ignite  an  exterior  wall.  All  of  these  mechanisms  depend 
essentially  on  direct  contact  (convection  and/or  conduction)  for  the  heating  and  ignition 
of  the  combustible  exterior  materials  of  the  structure.  Alternately,  hot  brands  from 
burning  trees  remote  from  the  property  can  be  carried  by  the  wind  and  blown  up-against 
the  exterior  or  even  into  homes  through  small  openings  to  ignite  the  structure.  Or,  finally, 
for  large  fires,  thermal  radiation  emitted  from  the  flames  can  be  transported  in  line-of- 
sight  directions  across  open  spaces  to  be  absorbed  by  and  ignite  the  structure. 

Of  these  possible  modes  of  fire  spread,  the  most  quantitative  understanding  exists  for  fire 
spread  by  ignition  from  thermal  radiation  [9,  10,  43-47],  For  example,  radiative  heating 
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of  a surface  is  simply  additive;  therefore,  the  radiative  heat  flux  at  a position  on  the 
surface,  from  two  identical  non-interacting  burning  objects  oriented  similarly  with 
respect  to  that  position,  is  just  tw  ice  that  of  one  of  them.  In  addition,  radiation  becomes  a 
more  dominant  transport  mechanism  as  fires  become  larger.  Crown  fires,  which  burn 
intensely  and  spread  rapidly  relative  to  other  wildland  fires,  are  well-known  examples  of 
fires  where  radiation  is  a dominant  mode  for  fire  spread.  The  fact  that  the  heat  transfer  by 
radiation  from  burning  trees  to  a structure  could  be  quantified  through  knowledge  of  the 
geometric  relationship  between  the  tw  o fuel  elements,  lead  to  the  selection  of  this  mode 
of  ignition  as  to  be  the  first  to  be  quantified  in  EcoSmart  [48-5 1 ].  This  was  done  with 
full  knowledge  that  more  often  structures  are  ignited  by  other  means  that  are  not  well 
quantified  as  discussed  above.  A simple  algorithm  was  constructed  to  determine  the 
thermal  radiation  transport  from  trees  to  structures. 

The  fire  spread  algorithm  can  be  used  to  examine  many  scenarios  for  fire  spread  from 
vegetation  to  structures.  The  algorithm  provides  means  to  calculate  the  radiative  heat 
exposure  to  structures  from  burning  trees  with  flames  located  at  arbitrary  positions  in 
3-dimensions  relative  to  structural  elements. 

Although  it  is  well  known  that  piloted  radiative  ignition  of  combustible  materials  depends 
on  both  the  magnitude  and  the  duration  of  the  incident  radiative  heat  flux  [48-51],  the 
model  described  here  introduces  only  a simple  critical  flux  assumption  for  ignition.  To 
be  consistent  with  other  capabilities  of  the  initial  version  of  EcoSmart,  the  burning 
duration  was  eliminated  from  the  algorithm  for  the  ignition  model.  Fire  spread  from 
vegetation  to  the  structure  is  assumed  to  occur,  therefore,  if  the  value  of  the  radiative  flux 
exceeds  a critical  value  for  ignition.  In  Section  5.1.6,  there  is  additional  discussion  of  this 
assumption,  and  a value  of  this  critical  flux. 

The  computational  relationships  that  form  the  basis  for  the  initial  fire  spread  model 
algorithm  (version  1 .0)  for  EcoSmart  are  summarized  below. 

5.1  SIMPLE  ALGORITHM  FOR  STRUCTURE  IGNITION  BY  BURNING 
LANDSCAPE  TREES 

The  algorithm  consists  of  six  steps:  1)  determine  the  burnable  mass  of  the  tree,  2) 
determine  the  flame  height,  3)  determine  the  peak  HRR,  4)  determine  the  fraction  of  the 
peak  HRR  emitted  as  radiative  energy,  5)  determine  the  fraction  of  the  radiative  flux 
emitted  by  the  flame  that  impacts  any  part  of  the  exterior  surface,  and  6)  determine  if  the 
criteria  for  ignition  of  the  exterior  surface  material  is  satisfied.  The  algorithm  has  been 
structured  to  utilize  the  burnable  mass  of  the  tree  as  a single  parameter,  to  which  burning 
properties  are  related.  Future  research  may  show  that  this  single  parameter  can  be  used  to 
generalize  burning  results  betw  een  many  different  types  of  trees.  As  this  study  examined 
only  one  type  of  tree,  however,  experiments  to  determine  additional  data  required  for 
generalization  of  the  results  to  other  tree  species,  w ithout  generation  of  additional  tree- 
bum  experiments,  is  not  recommended. 
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The  measured  mass  loss  during  burning  is  used  as  the  closest  approximation  to  the 
burnable  mass  available  for  the  trees  burned  in  this  study.  The  ratio  of  burnable  mass  for 
the  trees  to  the  pre-burn  mass  ranged  from  0.27  to  0.65. 

5.1.1  Burnable  Mass  of  the  Tree 

The  burnable  mass  of  the  tree  is  an  input  to  the  algorithm.  This  is  often  referred  to  as  the 
foliar  mass(m,  //  ;/.)  of  the  tree.  This  information  is  obtained  from  other  sources  within 

the  EcoSmart  software. 

5.1.2  Flame  Height 

The  flame  height  at  peak  burning  rate  is  related  to  the  crow  n height  of  the  tree  as  shown 
in  Figure  5.  This  data  is  complicated  by  the  fact  that  both  the  size  and  the  moisture 
content  of  the  trees  varied.  More  data  is  needed  to  understand  the  best  way  to  relate 
flame  height  to  the  tree  properties.  For  the  algorithm,  a fixed  ratio  of  flame  height  to 
crown  height  equal  to  2.0  is  used.  This  choice  is  supported  by  the  limited  data. 

Therefore: 

H flame  = 2.0  Hcr0wn  (2) 

where, 

H flame  [m],  length  of  the  flame  above  the  bole  of  the  tree, 

Hcrown  [m],  vertical  distance  from  the  bole  to  the  highest  branch. 

5.1.3  Peak  Heat  Release  Rate 

The  peak  HRR  is  well  know  n to  be  strongly  dependent  on  the  moisture  content  of  the 
burnable  mass.  In  a study  of  2 m Douglas-fir  tree  burning  conducted  by  Babraukas  [39] 
the  ratio  of  peak  HRR  decreases  exponentially  with  increasing  moisture  content  of  the 
needles  over  the  range  in  which  the  tree  could  be  ignited.  The  correlation  developed  by 
Babraukas  relates  the  peak  HRR  of  the  burning  tree  to  weight  of  the  tree  at  ignition.  For 
an  algorithm  to  be  used  in  EcoSmart,  the  data  must  be  presented  in  terms  of  burnable 
mass.  In  the  Douglas-fir  trees  burns  conducted  in  this  study,  the  pre-ignition  tree  weight 
was  approximately  twice  the  mass  loss  during  burning  (see  Table  3).  As  this  fraction  of 
the  initial  tree  mass  lost  during  burning  is  the  best  approximation  to  the  burnable  mass 
( yyi /l)liur ) available  from  measurement  made  in  this  study,  the  correlation  published  by 

Babraukas  has  been  multiplied  by  the  factor  2 to  adjust  for  use  with  the  burnable  mass 
parameter.  The  peak  HRR  for  the  algorithm  is: 


peak 


— 2 m 


5.S4-0.017A/ 


foliar 


where 


(3) 
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qpeak  is  the  peak  heat  release  rate  [kW] 
mJhliar  is  the  burnable  mass  of  tree  [kg] 

M is  the  moisture  content  of  the  needles  [in  percent  of  the  needle  mass] 


5.1.4  Radiative  Fraction 

The  fraction  of  heat  emitted  from  a fire  as  thermal  radiation  can  vary  widely  (see 
Quintiere  [10]).  Here,  it  is  chosen  to  be  a mid-range  value  of  0.35  (pp.  58-59  of  reference 
[10]  and  HUD  Report  [52]). 


*7  rad  0.  J)  5 Cj  peaj. 


(4) 


where, 

q rul  [kW]  is  the  peak  emitted  thermal  radiation  from  the  flame 

5.1.5  Radiative  Heat  Flux 

The  heat  flux  at  any  part  of  the  exterior  surface  of  the  structure  is  determined  by  a 
calculated  view  factor.  The  view  factors  are  based  on  radiation  from  a vertical  right 
cylinder  body  representing  the  flame  orientated  in  3-dimensions  relative  to  the  structure. 
The  height  of  the  cylinder  is  equal  to  the  flame  height  (H  flame ),  with  the  bottom  of  the 
cylinder  position  at  the  location  of  the  bole  of  the  tree.  The  calculation  of  view  factors  is 
given  in  Appendix  A. 

5.1.6  Ignition  Criterion 

The  ignition  criterion  for  the  structure  is  based  on  a critical  radiant  flux  condition.  There 
is  no  dependence  on  exposure  time.  A single  value  of  3 1 .5  kW/nr  is  used.  This  value  is 
based  on  criteria  developed  by  the  Department  of  Housing  and  Urban  Development  for 
safe  spacing  of  housing  from  potentially  large  fires  [52,  53],  This  value  is  relatively  high 
compared  to  minimum  heat  flux  values  needed  to  ignite  combustible  exterior  materials. 
However,  without  any  time  dependency  for  the  ignition  model  in  the  algorithm  and  using 
only  the  peak  radiant  energy  as  exposure,  this  value  is  suitable  for  current  use. 

6.0  IMPLEMENTATION  OF  THE  ALGORITHM  IN  ECOSMART 

The  fire  model  algorithm  was  selected  to  be  physics-based  and  to  harmonize  as  much  as 
practical  with  the  other  modules  of  EcoSmart.  Many  choices  that  were  made  in  the 
development  of  the  algorithm  for  ignition  of  structures  by  burning  vegetation  were  made 
specifically  to  work  within  the  limitations  of  the  EcoSmart  software  design.  There  is  a 
greater  quantitative  understanding  of  fire  spread  from  vegetation  to  structures  than  is 
implemented  in  the  softw  are.  Full  and  detailed  models  of  fire  development  on 
landscaped  properties  have  been  performed  with  the  NIST  developed  Fire  Dynamics 


19 


Simulator  (FDS)  software  [35-37],  Use  of  FDS  directly  in  EcoSmart  at  this  time  is  not 
practical  because  both  the  computational  resources  and  the  data  requirements  are  too 
extensive.  Using  the  algorithm  developed  where  ignition  of  structures  is  based  solely  on 
fire  radiant  exposure  should  only  be  regarded  as  one  indicator  of  fire  spread  potential 
from  burning  vegetation  to  structure. 

As  already  discussed,  the  initial  structure-ignition  algorithm  for  EcoSmart  does  not 
account  for  many  well  known  fire  spread  mechanisms  that  threaten  houses.  For  example, 
the  model  algorithm  does  not  contain  any  effect  of  wind.  It  is  well  known  that  fire  spread 
and  spread  rate  are  enhanced  by  wind,  see  for  example  Chapter  5 of  Quintiere  [10]. 
Never-the-less  as  w ind  effects  are  not  addressed  in  the  first  version  of  EcoSmart,  they 
are  not  included  in  the  structure-fire  algorithm.  Although  many  species  of  trees  and 
shrubs  will  be  available  for  use  in  the  EcoSmart  software,  it  was  found  that  very  little 
data  is  available  that  describes  the  burning  of  individual  trees  and  shrubs  in  a form  that  is 
needed  for  the  algorithm.  As  a result,  new  data  was  generated  by  NIST  for  burning 
Douglas-fir  trees.  As  the  vegetation  burning  portion  of  the  algorithm  was  developed 
solely  on  data  from  Douglas-fir  tree  bums  use  of  the  Version  1 .0  algorithm  with  other 
trees  is  not  recommended.  Additional  research  and  data  collection  is  needed  to  extend 
the  algorithm  to  other  tree  species. 

As  many  potentially  important  fire  spread  mechanisms  have  been  omitted,  conservative 
choices  that  favored  structure  ignition  have  been  made  for  the  fire  scenarios.  This  is 
customary  practice  in  safety  analysis  where  considerable  uncertainties  exist.  In  the 
algorithm: 

• All  vegetation  burns  at  the  maximum  burning  rate. 

• All  surfaces  ignite  if  exposed  to  a critical  flux  for  ignition. 

• Duration  of  burning  for  the  sources  of  heat  are  not  considered 

• Duration  of  exposure  for  the  materials  being  ignited  are  not  considered 


Since  this  version  of  EcoSmart  has  only  implemented  a model  of  structure  ignition  by 
thermal  radiation,  other  important  ignition  threats  to  structures  have  been  deferred  for 
later  development  and  implementation.  For  example,  the  initial  version  of  the  algorithm 
does  not  include  the  effects  of  ambient  winds,  burning  brands  and  heat  fluxes  from 
nearby  burning  structures,  for  example.  The  determination  of  the  relative  importance  of 
each  of  these  effects  is  scenario  dependent,  and,  therefore,  inclusion  of  the  effects  will 
require  additional  research,  data,  and  model  development. 

7.0  MODELING  COMPLEX  PHENOMENA 

A mathematical  model  of  any  phenomenon  can  be  regarded  as  a caricature  of  that 
phenomenon.  By  its  very  nature,  modeling  is  a creative  venture.  Just  as  a caricature 
depends  on  the  artist  for  its  rendition,  so  mathematical  modeling  depends  on  the  modeler 
for  its  mathematical  description.  If  a model  is  well  crafted  and  efficient,  it  will  capture 
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the  important  features  of  the  phenomenon  while  ignoring  less  relevant  ones  and  will 
produce  results  that  are  good  representations  of  the  phenomenon. 

Our  ability  to  construct  effective  mathematical  models  depends  on  our  understanding  of 
the  phenomenon.  If,  for  example,  we  do  not  understand  a phenomenon  conceptually,  we 
will  not  be  able  to  develop  a mathematical  representation  of  it.  As  an  aside,  however, 
any  attempt  to  construct  a mathematical  model  often  leads  to  a much  better  understanding 
of  research  needs  by  focusing  attention  on  developing  relationships  between  variables 
and  data  some  of  which  may  not  be  available. 

Mathematical  models,  even  of  the  same  phenomenon,  can  vary  widely  in  complexity,  as 
measured  by  several  criteria:  by  the  initial  research  and  development  effort  required  to 
construct  the  model,  by  the  computational  resources  needed  to  use  the  model  and  by  the 
data  required  for  input  and  validation  of  the  model.  Using  cost  as  the  distinguishing 
feature,  Figure  9 illustrates  different  classes  of  models  by  three  caricatures  of  a house, 
ranging  from  most  detailed  on  the  left  to  the  simplest  on  the  right.  We  refer  to  the  sketch 
on  the  left  as  a one-dollar  model,  the  sketch  on  the  right  as  a five  cent  model,  and  the 
center  sketch  as  a 25  cent  model.  The  price  of  each  model  is  used  here  to  convey  the  idea 
that  greater  model  complexity  generally  leads  to  higher  research  and  development  costs, 
more  extensive  computer-resource  requirements  and  greater  data  needs  both  for  input  and 
validation  of  the  model.  It  also  implies  that  greater  model  complexity  usually  yields  a 
more  faithful  and  easily  understood  representation  of  the  phenomenon. 


Complexity  is  neither  a reliable  gauge  of  model  accuracy  nor  effectiveness.  The  goal  of 
any  mathematical  model  is  to  produce  the  simplest  representation  of  that  phenomenon 
consistent  with  user  needs  to  understand  and  predict  its  behavior.  It  is  generally  true  that 
representing  a complex  phenomenon,  like  fire,  requires  a complex  model  (a  dollar 
model).  For  example,  a simple  (five  cent)  model,  while  often  adequate  to  illustrate  may 
be  inadequate  or  misleading  for  evaluations  of  details. 
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Figure  9 Three  caricatures  of  a house  representing  different  levels  of  complexity  for 
mathematical  models. 
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The  mathematical  model  algorithm  for  ignition  of  structures  used  in  version  1.0  of 
Ecosmart  is  distinctly  a 5-cent  class  model.  As  noted  earlier,  many  choices  were  made 
throughout  the  development  of  the  algorithm  to  simplify  the  model  in  order  to 
approximate  the  ignition  potential  of  a structure  due  to  nearby  burning  vegetation,  while 
requiring  minimal  computer  resources  and  minimal  input/validation  data.  During  a 
review  of  a pre-release  version  of  the  fire  portion  of  the  EcoSmart  software,  many  of  the 
model  limitations  were  identified  by  the  reviewers.  These  limitations  have  served  as  the 
starting  point  for  defining  research  needs  to  improve  future  versions  of  the  software. 

8.0  RESEARCH  NEEDS  TO  ADVANCE  THE  DEVELOPMENT  OF  PHYSICS- 
BASED  MODELS  FOR  WUI  AND  RESIDENTIAL  HOUSING  COMMUNITY 
FIRES 

As  discussed  above,  operational  models  are  regularly  used  in  practice  by  land-managers, 
and  fire-prediction  specialists,  for  example,  to  estimate  the  spread  and  spread  rate  of 
wildland  fires.  However,  these  models  are  limited,  as  noted  above,  and  badly  need 
updating  to  take  into  account  of  current  mathematical  modeling  methodology  and  the 
revolutionary  advances  in  computer  technology  over  the  past  few  decades.  Since  the 
wildland  fuels  simulated  in  the  operational  models  have  no  vertical  structure,  the  fire 
growth  and  spread  are  not  governed  by  the  conservation  laws  of  physics.  Therefore,  there 
is  no  hope  of  generalizing  these  models  to  address  the  problem  of  fire  spread  to  buildings, 
transitions  from  ground  to  crown  fires,  or  fire  spread  from  individual  tree  to  tree.  In  other 
words,  these  models  will  never  be  able  to  model  fires  at  the  wildland-urban  interface! 

When,  in  addition,  the  built  environment  becomes  involved  in  a fire,  as  in  the  Oakland 
and  Berkeley  Hills  fire  of  October  21,  1991,  or  more  recently  the  Los  Alamos  fires  in 
May  2000,  Summerhaven,  AZ,  fires  in  June  2003,  and  San  Diego,  CA  fires  in  October- 
November,  2003  these  models  are  ineffective.  The  operational  models  cannot  predict  the 
spread  of  fire  because  the  fuel  load  from  buildings  needs  to  be  treated  as  discrete  and 
three-dimensional.  Also,  in  many  cases  the  fuel  load  from  structures  is  much  larger  than 
the  surrounding  vegetation.  In  order  to  model  WUI  and  housing  community  fires, 
buildings,  as  well  as  large  individual  trees,  must  be  regarded  as  discrete  fuel  elements  on 
the  length  scales  of  interest.  At  a fundamental  level,  the  empirical  correlations  upon 
which  the  wildland-fire  operational  models  have  been  developed  are  no  longer  valid. 
Predictive  methods  are  needed  to  address  the  potential  fate  of  housing  communities  and 
individual  high  valued  structures  threatened  by  fire.  No  validated  predictive  models  of 
fires  in  an  urban  or  urban/wildland  setting  exist.  Research  is  needed  to  develop  three- 
dimensional  models  of  fire  phenomena  in  WUI  and  housing  communities.  Measurement 
of  three-dimensional  burning  characteristics  of  individual  vegetation  and  man-made 
structure  fuel  elements  are  needed  to  support  the  model  development. 

9.0  NIST  FIRE  DYNAMIC  SIMULATOR  AND  SMOKEVIEW  MODELS 

Over  the  past  25  years,  the  Building  and  Fire  Research  Laboratory  (BFRL)  at  the 
National  Institute  of  Standards  and  Technology  (NIST)  has  developed  a physics-based 
mathematical  and  computational  model,  known  as  the  Fire  Dynamics  Simulator  (FDS),  to 
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predict  fire  phenomena  inside  structures.  Over  time,  this  model  has  been  used  to  predict 
outdoor  buoyant  plume  behavior,  outdoor  fire  spread  from  structure  to  structure,  and, 
most  recently,  fire  spread  in  both  continuous  and  discrete  wildland  fuels.  FDS,  and  a 
NIST-developed  companion  code  for  visualization  called  Smokeview  have  been  well 
documented,  validated  using  data  from  large-scale  indoor  and  outdoor  experiments,  and 
are  used  by  hundreds  of  fire  protection  engineers  around  the  world.  The  model  equations 
and  their  computational  implementation  in  FDS  have  been  carefully  selected  to  be 
computationally  efficient.  Therefore,  while  FDS  is  both  computationally  and  data 
intensive,  simulations  requiring  on  the  order  of  a million  or  more  computational  cells 
generally  require  only  a current  high-end  PC  running  overnight.  Both  codes  and 
documentation  can  be  downloaded  free  of  cost  from  the  URL:  http://fire.nist.gov. 

The  generation  of  mathematical  models  represented  by  FDS  is  based  on  the  methodology 
of  Computational  Fluid  Dynamics  (CFD),  and  takes  advantage  of  the  revolutionary 
advances  in  computer  technology  with  correspondingly  reduced  costs  that  make  effective 
fire  modeling  capabilities  widely  available.  There  is  a parallel  between  the  current 
operational  models  for  wildland  fires  and  the  previous  generation  of  fire  models 
developed  by  NIST  in  the  1980’s.  These  “zone-models”  were  effective  in  computing  fire 
phenomena  in  structures,  but  only  as  average  conditions  on  a room-by-room  basis.  The 
need  for  engineers  to  resolve  greater  details  within  rooms  forced  a transition  to  the 
present  generation  of  FDS  fire  modeling.  To  advance  fire  modeling  of  WUI  and 
residential  communities,  a new  generation  fire  modeling  tools  must  be  developed  for 
general  use  by  land  managers,  community  authorities,  and  property  owners. 

At  this  time,  there  seems  to  be  little  to  no  research  activity  within  the  wildland-fire 
community  directed  at  developing  the  next  generation  of  three-dimensional  fire  models. 
NIST  has  begun  to  expand  the  capabilities  of  the  FDS  fire  model  and  Smokeview  fire 
visualization  software  to  predict  important  fire  phenomena  occurring  in  the  WUI  and 
residential  housing  communities.  This  is  a giant  modeling  undertaking.  To  estimate  the 
order  of  magnitude  of  the  research  task  to  develop  and  validate  a next- generation  fire 
model,  it  has  been  estimated  that  NIST  and  its  other  agency  partners  have  invested  about 
100  scientific  staff-years  in  the  development  of  the  present  FDS  and  Smokeview'  tools! 
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